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We show that the onset of laser multiple filamentation can be described as a critical phenomenon 
that we characterize both experimentally and numerically by measuring a set of seven critical ex¬ 
ponents. This phase transition deviates from any existing universality class, and offers a unique 
perspective of conducting two-dimensional experiments of statistical physics at a human scale. 
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The extensive study of phase transitions has triggered 
off substantial experimental and theoretical achieve¬ 
ments. Critical phenomena can be sorted out into uni¬ 
versality classes, each one being characterized by a set 
of critical exponents. However, due to the intrinsic ex¬ 
perimental difficulties arising when one needs to tweak 
control parameters, most of the theoretical predictions 
have only been tested numerically, leaving on^ a few ex¬ 
perimental verifications of generic models [l|, [2[ available 
to date. 

Population inversion at the core of laser phenom¬ 
ena has been the first example of an out-of-equilibrium 
second-order phase transition Q . Since then, as the avail¬ 
able optical power gradually increased throughout the 
years, laser physics opened the way to the field of non¬ 
linear optics, bringmg upfront the stucW of solitons Q, 
optical vortices [5|, l6| , or filamentation [7H9| . 

The latter phenomenon describes a self-guided propa¬ 
gation regime typical of high power lasers, and is initiated 
when the power exceeds a critical value Per, allowing the 
Kerr self-focusing phenomenon to overcome diffraction. 
The self-focusing laser pulse then reaches intensities al¬ 
lowing medium ionization. The resulting plasma as well 
as other higher-order effects in turn contribute to the 
laser defocusing, therefore initiating a dynamic balance 
between these two competing effects. 

When the input power largely exceeds Per, the beam 
breaks up into multiple cells that self-focus individually. 
This process is driven by optical turbulence [Tol-llij. in 
which transverse modulational instability seeded by the 
irregularities of the laser intensity profile as well as the in¬ 
homogeneities of the atmospheric refractive index results 
in the local nucleation of filaments across the beam pro¬ 
file. The multiple filamentation pattern steadily builds 
up as the pulse propagates. Progressively, each individ¬ 
ual filament reaches its clamping intensity by absorbing 
energy from the surrounding photon bath, that gets de¬ 
pleted. The beam therefore evolves from an initially noisy 
but homogeneous, fully-connected fiuence state towards 
a clustered transverse section, with isolated plasma fil¬ 
aments surrounded by independent islands of moderate 
fiuence. 


This process results in typical patterns featuring as 
much as 1000 filaments, that grow at the expense of 
the photon bath intensity. Their spatial concentration is 
limited by lateral interactions, which prevent them from 
packing closer to one another beyond a density of several 
cm-2 

Based on such experimental multiple filamentation 
patterns as well as numerical simulations, we show that 
during the beam propagation, a sharp transition oc¬ 
curs as the emergence of multiple filaments structures 
the beam into isolated high-fiuence clusters separated 
by moderate-fiuence areas. This evolution constitutes a 
phase transition, that we characterize by a set of seven 
critical exponents. Their values verify a hyper scaling uni¬ 
versal relationship predicted by renormalization group 
arguments thus confirming the existence of a crit¬ 
ical phenomenon. However, in spite of the similarity 
with percolation patterns, the critical exponents devi¬ 
ate from those of the corresponding universality class, as 
well as from any other previously characterized system. 
The multiple filamentation onset therefore exhibits a new 
and experimentally accessible universality class, exempt 
from the usual drawbacks of sample preparation or is¬ 
sues on the control parameters. Our findings may also 
pave the way towards an original theoretical framework 
for investigating the effect of the geometry of multiple 
filamentation on applications like laser-induced conden¬ 
sation 0 or the triggering of lightning [l8|, where the 
mechanism at a microscopic scale is strongly impacted 
by connectivity within the laser beam. 

Experimental fiuence patterns from a 100 TW 
beam recorded on photographic paper over 15 m of prop¬ 
agation (Eigure[T^-e) illustrate the formation of a multi¬ 
ple filamentation pattern. Erom a single connected clus¬ 
ter (panel a), the depletion of the photon bath splits the 
beam into several distinct clusters (panels c-e), colored 
by size in EigurelH Details of the image processing can 
be found in the Supplementary Material. 

A closer look at the cluster structure (Eigure [2^-c) 
shows that, during the propagation, the “strings” con¬ 
necting the filaments vanish, while small-size (~ 2 mm^) 
clusters survive around individual plasma filaments. 
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FIG. 1. Formation of a structured disconnected filamentation pattern from an initially smooth one in a 100 TW beam. 
Experimental record (a-e) and numerical simulation (f-j) evolution of the fluence profile, after propagating in air for (a) 1.25 m; 
(b) 2.5 m; (c) 5 m; (d) 10 m, and (e) 15 m. Numerical simulations are shown here for (f) 1 m; (g) 8 m; (h) 10.6 m; (i) 11.7 m 
and (j) 14.4 m. The fluence has been thresholded at the same value in all pictures, corresponding to 0.8TW/cm^ for numerical 
simulations. Clusters are colored according to their relative size in each image. 



FIG. 2. Close-up on the (a-c) experimental and (d-f) nu¬ 
merical multiple filamentation patterns of Figure [U at prop¬ 
agation distances (a) 5 m; (b) 10 m, and (c) 15 m for the 
experimental data, and (d) 10.6 m; (e) 11.7 m and (f) 14.4 m 
for their numerical counterparts. Clusters corresponding to 
laser filaments survive all along the propagation, even when 
the filaments begin to vanish at long distances. 


The transition from a single cluster connecting the 
opposite sides of the beam profile, to a disconnected, 
multi-clustered pattern, appears to occur sharply be¬ 
tween z = 5 and 10 m, reminding us of a phase transition 
in a percolation experiment. 

In order to overcome the intrinsically limited amount 
of experimental data and their finite size, we sampled 
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FIG. 3. (a) Experimental mean cluster size (surface) as a 

function of the lattice filling factor p. (b) Maximum cluster 
size with respect to the system size L. We disregarded the 
percolating cluster appearing when p exceeds pc = 0.5 from 
the surface calculation. 


square lattices of different sizes from each image using 
a coarse-graining procedure. In each map, we measured 
the filling factor p, defined as the ratio between the num¬ 
ber of occupied sites and their total number, the spatial 
correlation length as well as the mean cluster surface S. 
In the latter measurement, following the usual practice 
of percolation analysis Hi , the percolating cluster is ex¬ 
cluded. 

As displayed in Figure [3^1, the mean cluster size ex¬ 
hibits a maximum value when p reaches close to a criti¬ 
cal value Pc ~ 0.5. Close to such criticality, for spatially 
infinite systems, the phase transition theory predicts a 
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divergence of the spatial correlation length ^ as well as 
the mean cluster size S', defining two critical exponents 
V and 7 : 


^ |p-Pc| (1) 

S^\p-p,\-\ ( 2 ) 

In a finite system, ^ is bounded to the width L of the 
sampled squares. Near criticality, it can therefore be 
assumed to lie close to this value, so that the previous 
scalings rewrite as: 






\p-Pc\ 


( 3 ) 

( 4 ) 


Based on this assumption, we fitted the ^ jv — 1.79 ±0.1 
ratio from our samples (Figure [ 8 ) 3 ). This value is inde¬ 
pendent within ±0.1 from the parameters of the coarse- 
graining procedure. 

Our result displays multiple filamentation as a very 
promising practical experimental system to investigate 
phase transitions. Indeed, as the propagation distance 
2 : increases, the filaments progressively emerge from the 
noise, and the filling factor p monotonically and contin¬ 
uously decreases as the intensity profile structures itself. 
In addition to being directly observable, this continuous 
and monotonic decrease also provides an easily and finely 
tunable control parameter for the said phase transition. 
Such property contrasts with many systems, where ex¬ 
perimental constraints drastically limit the fine-tuning of 
the control parameters, and in many cases only allow 
measurement in a few conditions 2 qI 42^ at the cost of 
complex and costly sample preparation [l|. 

Willing to fully characterize the phase transition at 
stake here, we turned to numerical simulations in or¬ 
der to supplement the experimental observations with 
a large amount of realizations and continuous measure¬ 
ments along the propagation distance. 


We simulated the propagation of a flat-top beam, with 
1 TW/cm^ incident intensity. As the experimental data 
exhibits a high shot-to-shot reproducibility, we infer that 
the initial phase mask is mainly governed by the laser 
profile rather than the atmospheric turbulence. There¬ 
fore, we modeled this intrinsic perturbation by adding a 
10 % white noise on amplitude and a 0.2 radian amplitude 
noise on the complex phase onto the input profile, rather 
than using turbulence masks US- Both of these noises 
were made continuously differentiable at the numerical 
grid precision thanks to an interpolation method. 

Comparing panels a-e and f-j of Figure [H as well as 
close-up panels a-c and d-f of Figure[2l respectively, shows 
that the model results exhibit major similarities with the 
experimental patterns, even if the detailed arrangement 
of the individual filaments differs to some extent. In 
particular, the emergence of connected filament clusters 
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FIG. 4. Evolution of (a) the percolation probability, (b) the 
filling factor p, and (c) the correlation length as a function 
of the propagation distance z, defining the phase transition 
between the percolating and non-percolating states. 


is adequately reproduced, as well as their size distribu¬ 
tion, its evolution during propagation, and the subse¬ 
quent splitting of the initial clusters into smaller ones of 
2 mm^ around the location of each filament. 

The simulated patterns were processed the same way 
the experimental ones were. Figure 0] summarizes the re¬ 
sults of this analysis, averaged over 92 independent noise 
realizations. As z increases, the filling factor monoton¬ 
ically decreases, causing the breakdown the correlation 
length as well as the fragmentation of percolating fluence 
patterns into multiple sets of non-percolating variously 
sized clusters. 

Besides the scaling laws on the average cluster size and 
correlation length (Equations ([T]) and (ED), two supple¬ 
mentary ones are accessible (see details in Supplemen¬ 
tary Material). They concern the cluster strength Pqo 
( i.e., the probability that a site belongs to the percolat¬ 
ing cluster), as well as the the typical maximum cluster 
size Smax- As expected for a critical phenomenon: 

-Poo ~ (p-Pc)^, P>Pc (5) 

^ \P - Pc\~'^■ (6) 


As summarized in Table [H the more detailed statistical 
analysis brought by the numerical simulations allows one 
to measure these critical exponents. Three more expo¬ 
nents were deduced from scaling relations, as detailed in 
the Supplementary Material. The measurement of seven 
critical exponents in a single experiment is exceptional 
since one often only has access to a few ones [23|, |24|. 
Up to date only K. Takeuchi [H achieved such a compre¬ 
hensive characterization, for the specific case of directed 
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percolation. 

Critical Filament at ion Filament at ion 2D uncorrelated 


exponent 

experiment 

simulation 

percolation 

/5 

- 

0.13 ±0.02 

5/36 a 

! 0.14 

7 

- 

3.0 ±0.3 

43/18 

2.34 

V 

- 

1.6 ±0.1 

4/3 « 

1.33 

7/1/ 

1.79 ±0.1 

1.91 ±0.05 

43/24 s 

1.79 

<7 

- 

0.33 ± 0.02 

39/91 

^0.39 

a 

- 

1 

L) 

* 

± 

0 

to 

-2/3 w 

-0.67 

T] 

- 

0.09* ± 0.29 

5/24 a 

! 0.21 

S 

- 

24.9* ± 6.5 

91/5 = 

: 18.2 

Pc 

0.5 ±0.1 

0.48 ±0.01 

0.5927 


TABLE L Critical exponents measured from multiple filamen- 
t at ion experiments and simulations, and expected values for 
two-dimensional uncorrelated percolation. The asterisks indi¬ 
cate that a value has been derived from a hyperscaling rela¬ 
tion, and not from a direct finite-size scaling measurement. 

Surprisingly, in contrast with the visual intuition from 
the similarity between percolation and multiple filamen- 
tation patterns (Figured]), these values demonstrate that 
the optical phenomenon does not belong to the same uni¬ 
versality class as uncorrelated two-dimensional percola¬ 
tion. 

This discrepancy is illustrated at long distances by the 
survival of clusters of a typical size of 2 mm^ (also seen 
in the experiments), even when p decreases down to 0.3. 
These are clearly visible in blue in Figure [2j and appear 
as a “topological defect” in the cluster sizes distributions 
(Figure [5|). 

Physically, this survival of clusters of well-defined size 
stems from the robustness of plasma filaments once they 
are initiated. Indeed, the dynamic balance at the root 
of filamentation can, e.g. withstand turbulence strengths 
5 orders of magnitude beyond natural atmospheric val¬ 
ues [m. This feature clearly diverges from uncorrelated 
percolation expectations, where the cluster size distri¬ 
butions obey a scaling law np{s) = s~^g±[{p — Pc)s^] 
close to criticality. Obviously, this equation cannot fit 
the size distribution corresponding to multiple filamen¬ 
tation (Figure [5|). 

It therefore appears that the emergence of multiple fil¬ 
amentation constitutes a phase transition. The reliabil¬ 
ity of our results is assessed by checking the universal 
scaling relation z/d = 2/3 + 7 , where d = 2 is the geomet¬ 
rical dimensionality of our system. In our case, we find 
(2/3 + 7)/z/ = 2.0 ±0.3, ensuring that we indeed deal with 
a critical phenomenon. However, the set of seven criti¬ 
cal exponents that we were able to determine from the 
numerical and experimental data matches none of known 
universality classes, suggesting that the patterns of laser 
multiple filamentation defines a new class by itself. 

Let us underline the fact that the critical exponents 
determined here are fundamentally equilihrium quanti- 


10 '^ 10 “' 10 ^ 
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FIG. 5. (a) Evolution of the normalized cluster size distri¬ 
butions Us with respect to the propagation distance. A large 
percolating cluster is visible until approximately z = 10m. 
Beyond, a mode emerges around 2 mm^, corresponding to a 
robust size associated with individual filaments, (b) Normal¬ 
ized cluster size distribution for z = 11.7m with respect to the 
cluster area in log scale. The typical size corresponds to the 
spatial extension of the photon bath immediately surrounding 
a plasma filament. 

ties, therefore associated to a given equilihrium situation 
at each considered propagation distance. If the patterns 
were to be considered as analogues of a classical spin sys¬ 
tem in the microcanonical ensemble, each image would 
then be an equilibrium configuration corresponding to a 
given temperature, controlled here by the propagation 
distance. In spite of the time evolution of the electric 
field, the consistency of our results indicates that the evo¬ 
lution of the filamentation pattern is quasi-static. Such 
a steady evolution can be related to the aforementioned 
filament robustness. 

It is unexpected to witness the success of equilibrium 
finite-size scaling techniques and their non-trivial predic¬ 
tions on patterns extracted from the integration of a dif¬ 
ferential equation as well as from real experimental mate¬ 
rial. The reason is twofold: firstly, the observed patterns 
arise from a modulational instability, but nothing a pri¬ 
ori suggests that non-trivial scaling laws would actually 
apply. Secondly, the variation of the filling factor with 
the propagation distance is intimately related to a time 
evolution, contradicting the standard picture of equilib¬ 
rium statistical physics. 

In conclusion, we have shown, based on both exper¬ 
iments and numerical simulations, that the emergence 
of multiple filamentation in a high-intensity, ultrashort 
laser beam is a critical phenomenon. This opens the op¬ 
portunity to realize two-dimensional percolation-like ex¬ 
periments at a human scale and relate to future theoreti¬ 
cal developments. Moreover, the values of the associated 
seven critical exponents evidences that this phenomenon 
belongs to a previously unidentified universality class. 

Beyond providing a first prototype for this universality 
class, the propagation of ultrashort laser pulses also offers 
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a very rich experimental and numerical benchmark, since 
the control parameter, namely the filling factor p, can be 
continuously tuned and controlled from 0 to 1 via the 
propagation distance 

We also expect the recently proposed atmospheric ap¬ 
plications of filamentation to directly benefit from our 
findings. For instance, the triggering of lightning [Til 
requires a high conductivity along the laser beam and 
therefore relies on connectivity between the multiple fila¬ 
ments, and the phase transition evidenced in the present 
work may have crucial impacts on the efficiency of these 
processes. 

We acknowledge the Forschungszentrum Dresden- 
Rossendorf (FZD) contribution (groups of Pr. R. Sauer- 
brey and Pr. U. Schramm) to the experiments obtained 
with the DRACO facility. Fruitful discussions with 
R. Vasseur, H. Duminil-Copin and N. Berti are also grate¬ 
fully acknowledged. We acknowledge support from the 
European Research Council Advanced Grant “Filatmo”. 


NUMERICAL SIMULATIONS 

In order to achieve the numerical simulation of the 
multifilamentation phenomenon, we turned to a two- 
dimensional time-integrated model (2D+1) [ 2 ^. Using 
a factorization ansatz for the real electric field, namely 
£ = A{x^y^z) X x(t), one obtains an effective equation 
for the field envelope A\ 


propagation step 6r] proportional to IV^Imax order to 
resolve the sharp intensity peaks. We assessed that the 
resolution of the numerical grid ( 1.8 ym per pixel in the 
transverse direction) was sufficient, by checking the con¬ 
servation of the system hamiltonian when losses due to 
the ionization were neglected (i.e., setting u = 0). More 
specifically, we monitored the conservation of the hamil¬ 
tonian H: 

^ = iivx^ii^ - \ m\i + (13) 

where || • ||n = [/ | * denotes the standard L'^ norm. 

Energy conservation was indeed ensured within less than 
1 % error over the propagation distance. 

Initial conditions were chosen as close as possible to 
the experimental situation. The fiat-top beam was mod¬ 
elled by a radially symmetric fourth-order supergaussian, 
'0(r) = . The amplitude was calculated 

from the initial power Pin by 


'fpo 


7rcr^r(5/4) V 7 / 


(14) 


To ensure a continuously differentiable intensity and 
phase profile, the noise was first calculated on a coarse 
mesh, then interpolated on our 4 fold larger numerical 
grid. Intensity profiles are stored over the whole propa¬ 
gation length of 15 m for subsequent analysis. 


iafcon2|A|2-iy|^|2^-77|^|2K-2 

(7) 

where is the transverse Laplacian, is the central 
wave-vector, a and 7 parameters modeling the delayed 
medium response, 77.2 the Kerr non-linear refractive index, 
and the multiphoton absorption coefficient {K=S 

in air). After rescaling the transverse dimensions by 

l/\/ 2 /co {"f/{oLkQn 2 )^) and introducing the con¬ 


venient set of dimensionless variables 

■q = (8) 

i. = {akon2l^-^) , (9) 

%l> = A{akon2h)~~^^‘'‘^^~‘^\ ( 10 ) 

eq. © rewrites under a compact generic form: 

+ AxtA + IV’P'0 + P(|'0|^)'0 = 0, (11) 

where the underlying physics is described by the func¬ 
tion F 

F(|V;|2) = -|V;|2^ + ii.|V;|2^-2. (12) 


We performed the numerical integration of eq. ([TT]) us¬ 
ing a split-step Eourier method and set an adaptative 


RECORDING OF EXPERIMENTAL PATTERNS 

71, 

As described in detail earlier M, experiments were 
performed using a titanium:sapphire chirped-pulse am¬ 
plification providing 3 J, 100 TW pulses of 30 fs dura¬ 
tion, at a repetition rate of 10 Hz and central wavelength 
of 800 nm. The beam was launched collimated in air 
with a diameter of approximately 10 cm, through a 6 
mm thick fused silica window, the dispersion of which 
was pre-compensated by adequately adjusting the grat¬ 
ing compressor of the laser system. The multiple fila¬ 
mentation of the beam was characterized by single-shot 
burns on the back side of photosensitive paper (Kodak 
Linagraph 1895). The images were then scanned into a 
grayscale colormap and processed as described in detail 
below. 

IMAGE PROCESSING AND ANALYSIS 

The beam fiuence profiles, whether experimental or 
numerical, corresponding to a propagation distance z, 
was handled as a grayscale image. The fiuence was first 
thresholded to a fixed value slightly below the initial fiu¬ 
ence level of the initial beam (0.8 TW/cm^ in the nu¬ 
merical simulations). We then coarse-grained the result- 


d^A = — 
2ko 
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ing black and white images on square lattices of 16 x 16 
pixels. 

Each bin was declared active or inactive (respectively 
filled or empty) depending on the proportion within the 
bin, of pixels above the fluence threshold from the original 
profile. We checked that the results were independent 
from the cell dimension as long as the coarse-graining 
was relevant. 

In the case of the experimental data, where only one 
realization was available, a statistical analysis was made 
possible by excluding the regions where diffraction pat¬ 
terns are clearly visible at the bottom of the beam pro¬ 
file, and by randomly sampling square systems of differ¬ 
ent sizes over the spatial region of interest (10 x 10 to 
64 X 64). In the case of numerical results, the large num¬ 
ber of realizations obtained at each distance allowed us 
to focus, for each of them, on the central square region 
of approximately 6 cm x 6 cm in the 50 TW case and 
8 cm X 8 cm in the 100 TW case, corresponding to the 
beam waist, on which we conducted all of our numerical 
measurements. 

We define a cluster as a set of active sites connected 
through nearest-neighbor interaction. Following [l9|, we 
define percolation when a cluster spans the lattice in both 
the horizontal and the vertical directions, in order to en¬ 
sure that one percolating cluster can exist at most. We 
calculate the percolation probability n(z), at any dis¬ 
tance z, as the average number of runs which percolate 
at this distance, divided by the total number of runs. 
We also investigated the the maximum cluster size Smax, 
the cluster strength Pqo (defined as the ratio between the 
largest cluster size and the system surface L^), the fill¬ 
ing factor p, the average cluster size P, as well as the 
correlation length The three latter observables are re¬ 
spectively defined as: 


p = y^jns 

(15) 

S=1 




S=-ys^n, 

Ppl 

(16) 

^2 yr'^gir) 

^ Erdir) ’ 

(17) 


where g{r) is the probability that two active sites sepa¬ 
rated by a distance r belong to the same cluster, i.e., are 
continuously connected by active sites. Us is the clus¬ 
ter size distribution, normalized by the system size L. 
Note that these values can be calculated for any given 
propagation distance z, or, alternatively, for any filling 
factor p. 

In the latter case, one has to take into account that, 
for each realization of the pulse propagation, the values 
of p fluctuate for a given z. Therefore, in order to per¬ 
form consistent ensemble averages, we weighted the ob¬ 
servables by assuming that p was distributed around its 


average following a Gaussian law of variance 
at each step, so that if one has n samples for a given 2 :, 
the average of any observable O reads 

{0{z)) oc ^ 0 ^e-(P.-p(^))V( 2 <T=(^))^ ( 18 ) 

/C=l 

where the normalization factor has been omitted for the 
sake of clarity. 


FINITE-SIZE SCALING 


Given the finite size of our samples, we determined the 
critical exponents /3, 7 , and a governing the evolution 
of Poo, *5, and Smax, respectively, close to the phase 
transition, by using finite-size scaling techniques. Such 
techniques rely on the ansatz that close to the percolation 
threshold, the correlation length is of the same order of 
magnitude as the system size ^ ^ L. 

The first step is to determine the phase transition 
threshold Pc. Let us for that purpose define the per¬ 
colation cumulant [Tol as the fraction 11 of configurations 
featuring a percolating cluster. Since p —Pc ^ L~'^ close 
to the percolation threshold, the scaling hypothesis al¬ 
lows one to write 


n{pr~^Pc,L) 


L- {Pc-p) 


(19) 


where is a generic scaling function. Note that the def¬ 
inition of n is only valid on finite lattices. It converges 
towards a step function at p = Pc when L 00 . Plot¬ 
ting n with respect to p for different system sizes exhibits 
a unique crossing point (Figure [ 6 ^), therefore yielding 
the value of Pc = 0.48 ± 0.01. Plotting 11 with respect 
to (pc — p) allows one to obtain a first estimation of u 
by curve collapse. Figure [Ob shows the determination of 

z/= 1.6 ±0.2. 

The exponent [3 defines the cluster strength behaviour 
as Poo ^ {p ~ PcY foi* P ^ Pc- Within the finite-size 
scaling hypothesis, introducing the proper function poo, 
the cluster strength is given by: 


Poo{p ^ Pc^L) = L -goo L-{p-pc) 


( 20 ) 


so that plotting PooL- with respect to p (Figure [^) 
should give the same unique crossing point as the perco¬ 
lation cumulant. Indeed, we find Pjv = 0.08 ± 0.01,and 
Pc = 0.48 ± 0.01, remarkably consistent with the previ¬ 
ous determination. Plotting the rescaled order parameter 
with respect to P^/^(pc —p) allows one to refine the previ¬ 
ous determination of v = 1 . 6 ± 0.1 by seeking a curve col¬ 
lapse (Figure [7b), from which we deduce /3 = 0.08 ± 0 . 02 . 


Finally, we determined 7 relying on the average cluster 
size. Using the finite-size scaling ansatz, S is expected 
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FIG. 6. (a) Percolation cumulant 11 with respect to the 

filling factor p. The expected unique crossing point de¬ 
fines Pc = 0.48 ± 0.01. All the curves result from the av¬ 
erage of 92 independent runs, (b) Percolation cummulant It 
as a function of the rescaled variable {pc — p)L^^^^ yielding 
z/?^1.6±0.1by curve collapse. 



FIG. 7. (a) Rescaled cluster strength Poo dependence on the 
filling factor p. The ratio ^/v — 0.08 ± 0.01 yields the same 
value for pc as its previous estimation using It. (b) Rescaled 
cluster strength Poo as a function of the rescaled parame¬ 
ter (jpc — p)L^^^^ allowing one to reduce the uncertainty on 
the correlation exponent, z/ = 1.6 ± 0.1. 


to increase as P? close to pc- Given the value of we 
obtain 7 = 3.0 ± 0.3 (Figure [ 8 ]). 

Finally, the exponent cr characterizes the mass diver¬ 
gence of the largest non-percolating cluster when ap¬ 
proaching criticality, namely Smax ^ \p ~ Pc\~^• Finite- 
size scaling ansatz therefore implies that 


<^max ijP ^ Pcf -^) — Qyi L ^ {j)q p) 


( 21 ) 


where, again, ^fM is a scaling function. Using v = 1.6, we 
obtain a data collapse for cr = 0.33 ± 0.02, as displayed 
on Figure [H 

The estimation of the exponents cr, S and z/ can be 
done using the hyperscaling relations a = 2 — dv (Joseph- 
son’s identity), 7 = (2 — 7 ^)zr, and 5 — 1 = ^^ (Widom’s 
identity). They lead to cr = —1.2 ± 0.2, 5 = 24.9 ± 6.5 
and 77 = 0.09±0.27. Note that the value of cr is consistent 
with Rushbrooke’s hyper scaling relation cr± 2 / 3±7 = 2 . 


FIG. 8. (a) Average cluster size S with respect to p. Note 
that we excluded the percolating cluster from the calculation 
for p > Pc. (b) Average cluster size with respect to the system 
size L, allowing us to find 7/z/ = 1.91 ± 0.05. 



FIG. 9. Rescaled largest non-percolating cluster as a func¬ 
tion of the rescaled parameter {pc — p)L^^^. Seeking a curve 
collapse leads to the determination of cr = 0.33 ± 0.02. 
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